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ABSTRACT 

We have investigated the post-merger signatures of red-sequence galaxies in rich Abell clusters at 
z < 0.1: A119, A2670, A3330 and A389. Deep images in u' , g' , r' and medium-resolution galaxy 
spectra were taken using MOSAIC 2 CCD and Hydra MOS mounted on a Blanco 4-m telescope at 
CTIO. Post-merger features are identified by visual inspection based on asymmetric disturbed features, 
faint structures, discontinuous halo structures, rings and dust lanes. We found that ~ 25% of bright 
(M r < —20) cluster red-sequence galaxies show post-merger signatures in four clusters consistently. 
Most (~ 71%) of the featured galaxies were found to be bulge-dominated, and for the subsample 
of bulge-dominated red-sequence galaxies, the post-merger fraction rises to ~ 38%. We also found 
that roughly 4% of bulge-dominated red-sequence galaxies interact (on-going merger) . A total of 42% 
(38% post-merger, 4% on-going merger) of galaxies show merger-relate d features. Compared to a 
field galaxy study with a similar limiting magnitude ( van Dokkum|[2005[ ) , our cluster study presents 
a similar post-merger fraction but a markedly lower on-going merger traction. The merger fraction 
derived is surprisingly high for the high density of our clusters, where the fast internal motions of 
galaxies are thought to play a negative role in galaxy mergers. The fraction of post-merger and on- 
going merger galaxies can be explained as follows. Most of the post-merger galaxies may have carried 
over their merger features from their previous halo environment, whereas interacting galaxies interact 
in the current cluster in situ. According to our semi-analytic calculation, massive cluster haloes may 
very well have experienced tens of halo mergers over the last 4-5 Gyr; post-merger features last that 
long, allowing these features to be detected in our clusters today. The apparent lack of dependence 
of the merger fraction on the clustocentric distance is naturally explained this way. In this scenario, 
the galaxy morphology and properties can be properly interpreted only when the halo evolution 
characteristics are understood first. 

Subject headings: catalogs - galaxies: clusters: individual (Abell 119, Abell 2670, Abell 3330, Abell 
389) - galaxies: evolution - galaxies: formation - galaxies: structure 
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1. INTRODUCTION 

The formation of massive early-type galaxies in the 
universe is still in question. The hierarchical galaxy 
formation scenario is widely accepted at present, which 
predicts that massive galaxies form through hierarchical 
galaxy mergers. If each galaxy merger induces star for- 
mation and consequent stellar mass growth in the united 
galaxy, it would lead to a large scatter in the age and 
metallicity of stellar populations in massive early-type 
galaxies. However, the observational characteristics of 
early-type galaxies, such as their red colors tightly cor- 
related in optical color-magnitude relations (CMRs) and 
their high a— elements ratios, imply that most of the stel- 
lar contents in massive early-type galaxies formed in a 
short timescale at an early epoch [z > 1). Furthermore, 
the red-sequenc e in the optical CMR appears to be estab- 
lished by z ~ 1 ( Tanaka et al.|2005 1 , and recent observa- 
tions with 8-10m-class telescopes reveal the appearance 
of massive red-sequence galaxies (M > IO^Mq) at z ~ 2 
3 JKodama et al.|2007j|Drory fc Alvarez|2008| |Kang fc 
li. 



These observational results suggest that most 
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massive early-type galaxies had almost completed their 
star formation and mass aggregation by z ~ 1 and then 
virtually passively evolved. 

There are other observational clues which indicate a 
significant increase in the stellar mass density in massive 
red-sequence galaxies since z ~ 1, which is not allowed 
according t o the p assive evolution of blue galaxies alone. 



Bell et al. (2004) showed that the B-band luminosity 
density ot massive red-sequence galaxies does not evolve 
much in the redshift range of < z < 1.1. Because 
the i?-band light should dim as stars get old, they ar- 
gued that some young populations should be provided 
regularly from blue galaxies to the red-sequence during 
this period. However, the simple fading of massive blue 
galaxies cannot be applied, as the brightest red galaxy is 
always brighter than the brightest blue galaxy through- 
out the redshift range. After all, they suggested that 
galaxy mergers are an important process in the forma- 
tion of luminous re d galaxies sin c e z ~ 1. The result was 
supported again by Faber et al. ( 2007 1 . 

Residual star formation (RSF) in early-type galaxies 
provides another clue about galaxy formation. RSF 
has been extensively studied, recently using the UV 
data from space telescopes such as the Galaxy Evolu- 
tion Explore (GALEX) and the Hubble Space Telescope 
(HST) . The latest work on the UV upturn phenomenon 
of early-type galaxies showed that RSF fraction among 
cluster elliptical galaxies is as high as ~ 30%; it is even 
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higher in the field environment ( Yi et aL|2011 ). Several 
studies have claimed that the RSF in early-type galaxies 
is most likely related to galaxy merg ers or interactions 
( |Kaviraj et~aLl|2007| |Kaviraj||2010b| . Because the UV 
bright phase or a young stellar population lasts only ~ 1 
Gyr, we can estimate that the RSF detected at less than 
z = 0.1 is stimulated at a relatively low redshift, z ~ 0.2 
- 0.3. Again, this is indirect evidence of the substantial 
fr equency of ga l axy m ergers given a low redshift. 



van Dokkum 
around field 



(2005|) showed tidally disturbed features 

0.1 with opti- 



liptical galaxies at z 
cally deep images from NDWFS (NOAO Deep Wide 
Field Survey) and MUSYC (MUlti-wavelength Survey 
of Yale and Chile). That paper suggested that ~ 70% 
of field elliptical galaxies were assembled through dry 
mergers (i.e., mergers of gas -poor, bulge-dom inated sys 



terns) in the recent epoch. Kaviraj (2010a) also found 



disturbed features as well as dust lanes from 25% of lu- 
minous early-type galaxies (M r < —20.5, z < 0.05) in the 
'Strip 82' fields of Sloan Digital Sky Survey (SDSS) data, 
which have much longer integration times than general 
SDSS fields. These are very meaningful results because 
mass growth via continuous galaxy mergers were up to 
that point mostly supported observationally by statisti- 
cal analyses of photometric colors instead of by visual ev- 
idence of galaxy mergers. Moreover, the merger fraction 
is often estimated from gal axy pair fractions, considering 



them as pre -mergers (e.g., Bundy et al. 2009 de Ravel 



et al.|l2009D . Although van DokkumT "( |26o5p an 
raj| ( |2U10a| did not deal with complete volume-limited 
samples, they searched a substantial number of galaxies 
with unprecedentedly deep optical images and presented 
direct evidence of galaxy mergers related to the forma- 
tion of red, early-type galaxies. 

Our interest moved to galaxy clusters, which are dom- 
inated by massive early-type galaxies. It has been ex- 
pected that galaxy clusters are not a likely environment 
for frequent merger events to take place due to the high 
peculiar motions of the galaxies within the cluster. More- 
over, it is commonly held that massive galaxies in clus- 
ters formed faster at an earlier epoch than did galax- 
ies in a field environment (Gunn & Gott 1972; Dressier 
1980; Tanaka et al. 2005 and references in it). There- 
fore, it was presumed that it would be difficult to find 
post-merger signatures from massive early-type galaxies 
in galaxy clusters. Naturally, there has been less effort 
to find post-merger signatures among cluster galaxies. 

To find post-merger signatures of early-type galaxies 
in galaxy clusters, we carried out optical deep imaging 
and multi-object spectroscopic observations of four rich 
Abell clusters at 0.04 < z < 0.11 with the Blanco 4-m 
telescope at the Cerro Tololo Inter-American Observa- 
tory (CTIO). We picked four rich Abell clusters (A119, 
A2670, A3330 and A389) as our targets. They have been 
observed by GALEX at least in the medium-depth imag- 
ing mode (MIS: Medium Imaging Survey) and are opti- 
mally visible from the CTIO. Their GALEX NUV expo- 
sure times range from ~ 1 hours for A119 to ~ 17 hours 
for A3330. UV light is very sensitive to the existence of 
young stellar populations and thus plays an important 
role in the study of the recent star formation history 
that may have been caused by mergers. The UV prop- 
erties will be published in an upcoming paper following 
this one, in which we initially present the optical colors 
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Fig. 1.— The field-of- views of the MOSAIC2 and Hydra data 
are superimposed over the DSS image for each cluster. The square 
in each panel indicates the field of the stacked deep image while 
the gray-dashed circles represent the observed Hydra FOVs. Their 
consistent FOVs enabled us to perform an almost complete spec- 
troscopic survey of the cluster galaxies. 

and morphologies of cluster galaxies. We performed a vi- 
sual inspection of all the red-sequence galaxies brighter 
than M r i — —20 in each galaxy cluster using deep optical 
images. This result will enable us to compare the recent 
galaxy merger histories between cluster and field environ- 
ments. Morphological indexes such as the bulge-to-total 
ratio and asymmetry of the galaxies were also measured 
to examine their morphological characteristics. 

In this paper, the observations and data reductions are 
presented in Section [2] and Section [3j Section [I] describes 
the galaxy sample selection process and Section [5] shows 
the scheme and the result of the visual inspection using 
the deep optical images. Morphological examination will 
be presented in Section [6j We will conclude this paper 
with a summary and discussion in Section [7] 

2. OBSERVATIONS 

Deep optical images of the Abell clusters were taken 
with the MOSAIC 2 CCD mounted on the Blanco 4- 
m Telescope at CTIO. Medium-resolution spectra (R ~ 
1000) of galaxies in the cluster fields were acquired with 
Hydra, a multi-object spectrograph on the Blanco. Ob- 
servations were performed three times through 2006 - 
2008. The target clusters and acquired data are summa- 
rized in Table [I] Only the short exposure image (t exp — 
600 sec) of A389 was taken separately in early 2009 with 
the same telescope and instrument. MOSAIC 2 consists 
of eight 2Kx4K CCDs (therefore, the overall dimension 
is 8Kx8K pixels), and Hydra utilizes a 2Kx4K CCD. 
The pixel scale of MOSAIC 2 is 0.24"/pixel after the as- 
trometry is done. MOSAIC 2 has a 36' x 36' field-of-view 
(FOV) and Hydra has a circular FOV that is 40' in di- 
ameter. The FOVs of MOSAIC 2 and Hydra are plotted 
over Digitized Sky Survey images in Figure IT] for each 
cluster. The square in each panel indicates the region 
of the stacked deep image while the gray-dashed circles 
represent FOVs of Hydra observations. Their consistent 
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TABLE 1 

CTIO OBSERVATIONS SUMMARY 



Cluster ID 


RA(J2000) 


Dec(J2000) 


Instrument 


Filter 




# of pointing 13 


Date 


Abcll 119 


00:56:21.4 


-01:15:47 


MOSAIC 2 


u' 


6000 / 1200 


5 


27,29 Oct. 2008 


(2 ~ 0.044) 








<j 


5760 / 120 


8 












r' 


5760 / 120 


8 










Hydra 




9000 


4 


o— 11 iNov. zUUo 


Abell 2670 


23:54:10.1 


-10:24:18 


MOSAIC 2 


vl 


4800 / 700 


6 


13,15 Nov. 2007 


(z ~ 0.076) 








g' 


5100 / 200 


6 












r' 


6000 / 60 


6 










Hydra 




04UU 


Q 

o 




Abell 3330 


05:14:40.0 


-49:03:15 


MOSAIC 2 


u' 


4600 / 460 


10 


26 Nov. 2006 


(z ~ 0.089) 








g' 


5000 / 200 


10 












r' 


6150 / 100 


15 










Hydra 




14400 / 3600 c 


3/6 


27-29 Nov. 2007 












/ 1-3 Dec. 2006 


Abell 389 


02:51:31.0 


-24:56:05 


MOSAIC 2 


u' 


6000 / 1200 


5 


27,29 Oct. 2008 


(z ~ 0.113) 








<j 


7200 / 1200 


6 












r' 


7200 / 600 


6 










Hydra 




10800 


3 


6-11 Nov. 2008 



a The first numbers are the total exposure time of the stacked images and the second numbers represent the short 
exposure images that were analyzed. 

b The numbers indicate the number of dithering points for imaging and the number of fiber configurations for 
spectroscopy for each target field. 

c Spectroscopic redshifts of some A3330 galaxies were complemented by the spectra taken in 2006. 



FOVs enabled us to perform an almost complete spec- 
troscopic survey of the cluster galaxies effectively. 

All of the deep images were taken on dark nights (moon 
illumination < 0.4, the average is 0.17) while the spec- 
troscopic observations were performed on bright nights. 
The imaging nights were usually photometric and the av- 
erage seeing was less than 1" during the entire period of 
observation. We obtained the deep images with v! , g', r' 
filters in the Blanco telescope. In order to fill the gap be- 
tween the eight chips in the MOSAIC 2 CCD, more than 
five points of dithering was applied for each filter. The 
maximum number of dithering points was 15 for the r' 
band deep image of Abell 3330. The u'g'r'i 'z' standard 
stars in Chandra Deep Field South(CDF-S, Smith et al. 
2003 1 were also observed every night to standardize the 
photometric data. 

Hydra observations were performed with a KPGL3 
grating centered at 5205A. A blocking filter, GG385, was 
also used to suppress the bluer wavelength below 3850A. 
The resolving power was R ~ 1100 at 5500A and the dis- 
persion was 1.40A/pixel after binning two pixels along 
the dispersion axis. Calibration frames were taken with 
a "Penray" lamp for wavelength calibration and with a 
"Quartz" lamp as an aperture reference. We used the 138 
large fibers available in Hydra. The fiber width was 300 
jum, which correspond to 2" in the sky. In a fiber config- 
uration, ~ 100 fibers were assigned for galaxies while ~ 
10 fibers were allocated to blank sky for the sky spectra. 
Three to four fiber configurations were designed for each 
galaxy cluster for the possible cluster members not to 
be missed by fiber collisions and the limited number of 
fibers. Every fiber configuration was also taken through 
three separate exposures for median combine for data 
reduction purposes. 

3. DATA 

3.1. MOSAIC 2 photometric data 



3.1.1. Data reduction 

The MOSAIC 2 data are mainly reduced using the 
mscred package in IRAF. All of the calibration data de- 
scribed in this section were obtained for every run and 
were applied to the object frames observed during the 
same run. First, we updated the characteristic informa- 
tion of each CCD chip. New crosstalk files were gener- 
ated with the xtcoef f task using shallow object frames 
which have saturated stars upon a small sky background 
for the CCD pairs of MOSAIC 2. Doing this corrects the 
ghost images of saturated stars reflected onto each other 
between a pair of CCDs. Gains and readout noises of 
eight chips were computed with the mscf indgain task. 
In addition, the bad pixel mask of each chip was revised 
manually via visual inspection. We then updated all of 
the FITS files of the same observing run with the new 
information. An IRAF task cedproc handled all of the 
pre-processes at once, such as crosstalk correction, bad 
pixel mask correction, overscan strip correction, image 
trim, bias correction, and flat-fielding. The combined 
dome flats were used for the flat-fielding of g' , r' band 
images, while the combined sky flat was used for the u' 
band images. 

Once the pre-process was done, we found a smooth 
background gradient in the g' , r' band deep images. Il- 
lumination correction was carried out for each band using 
the dark-sky flats generated from the object-subtracted 
deep images. After illumination correction, the deep im- 
ages showed ±2% deviation in their background over the 
field of each frame. Cosmic rays were removed from all 
of the long-expos ure images (t eX v ^ 500 sec) using the 
lacosmic routine ( van Dokkum|2001 1 . Satellite trails oc- 
casionally appeared in the long exposures. They were 
erased using the satzap task. 

Object frames were stacked to a single deep image in 
each band. The stacking procedure followed the gen- 
eral guide of the National Op tical Astronomy Observa- 
tory MOSAIC data reduction ( jJannuzi, Claver fc Valdes 
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2003 1 . Astrometry was carried out very carefully for ev- 
ery chip using the msctpeak task in IRAF to achieve 
positional errors of less than 0.3 " in RMS (Readers are 
referred to Valdes 2000 for details). 

3.1.2. Standardization 

Standardization equations were derived for each ob- 
serving run with u' g'r'i' z' standard stars in CDF-S 
( Smit h et al.| 2003 ). Basically, optical photometry in the 
band had to be performed using shallow images 
100 sec), as we had intentionally saturated the 



9 , r 

galactic centers in the deep images to see the faint struc- 
tures around the galaxies. Only for A389, the farthest 
galaxy cluster among our target clusters, did we use a 
single deep image (t exp = 1200 sec) and a relatively short 
exposure (t exp = 600 sec) for g' and r' band photometry, 
respectively. The standardization equations were estab- 
lished for the nights in which the short exposures were 
taken. 

The eight chips of MOSAIC 2 CCD have slightly differ- 
ent characteristics in terms of photon efficiency There- 
fore, great care was required when measuring the bright- 
ness from different chips for consistency. For this rea- 
son, we tried to allocate standard stars to the greatest 
extent possible in one chip and derived the standard- 
ization equations using the standard stars in it. Chip 
6 of the MOSAIC 2 CCD was dedicated to this work 
because it has the deepest electron-capacity among the 
eight chips. Stars with small photon-statistical errors (< 
0.03 mag) were used in the calculations. As a result, 7 
~ 9 stars were available for the standardization of the 
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band data and 6 stars were available for the 



band data a mong the 22 standard stars listed in Smith 
et al. ( 2003 1 . Instrumental magnitudes were measured 



with a fixed aperture dia meter, 14'/86, which is adopted 



from Smith et al. (2003). Because the aperture diame- 
ter is already larger than 7xFWHM in our data (typical 
FWHM < 1"), we did not consider aperture correction 
in the photometry of standard stars. 

Standardization equations were established as follows: 

(1) 

(2) 

(3) 



u ' = u inst + a w + b u >X u , + c u ,(u' - g') 
9* = 9i„st + a g' + h 9' x g' + c 9'(9' - r') 



The u'g'r' magnitudes in the above equations are the 
standard magnitudes and the instrumental magnitudes 
are defined as mi nst — —2.5 x log(counts(ADU/sec)) + 
25. X is an effective airmass of each object frame. Ini- 
tially, the basic offsets between the standard magnitudes 
and the instrumental magnitudes were calculated and 
corrected for the instrumental magnitudes. The mag- 
nitudes were then fitted versus the airmass. The fitting 
function provided the relationship with which to deter- 
mine the airmass term in the standardization equations. 
In the same manner, we fitted the airmass-corrected mag- 
nitudes versus the optical colors and determined the color 
terms in the equations. Weighted-mean values are used 
throughout the fitting process. The coefficients of stan- 
dardization equations are presented in Table [21 

The standardization equations were appliedto the in- 
strumental magnitudes from short exposure images to 
determine the standard magnitudes of galaxies. We then 



TABLE 2 

Standardization coefficients 



Date 


Filter 


a 


b 


c 


2006 Nov 26 


vl 


-1.478±0.064 


-0.441±0.009 


0.126±0.048 




g' 


1.002±0.012 


-0.168±0.005 


0.002±0.020 




r' 


0.982±0.005 


-0.088±0.002 


-0.018±0.007 


2007 Nov 13 


it' 


-1.385±0.086 


-0.399±0.018 


0.078±0.064 




g' 


1.028±0.021 


-0.139±0.012 


-0.049±0.029 




r' 


1.001±0.012 


-0.087±0.007 


-0.032±0.014 


2008 Oct 29 


u' 


-1.288±0.064 


-0.526±0.011 


0.095±0.048 




g' 


1.119±0.013 


-0.228±0.007 


-0.113±0.020 




r' 


1.054±0.006 


-0.155±0.003 


-0.041±0.009 



TABLE 3 
The average of /c-correction terms 



Cluster 


K„/ 


K g , 


K r , 


A119 


-0.20±0.03 


-0.11±0.02 


-0.04±0.01 


A2670 


-0.29±0.06 


-0.19±0.03 


-0.07±0.02 


A3330 


-0.34±0.06 


-0.24±0.04 


-0.09±0.02 


A389 


-0.42±0.06 


-0.32±0.05 


-0.13±0.03 



derived the standardization constants for the u'g'r' deep 
images using the standard magnitudes determined ear- 
lier. The offsets between the standard magnitudes and 
the instrumental magnitudes from the deep images were 
calculated with non-saturated stars using aperture pho- 
tometry. A typical standard deviation of the standard- 
ization constants (the average of magnitude offsets) was 
< 0.001 mag. 

3.1.3. Photometry 

The optical magnitudes of the galaxies were measured 
with SExtractor from th e MOS AIC 2 data and calibrated 
as described in Section |3. 1 .2 
Auto-Magnitu de (MAG_AuT0 



Essentially, w e adopted 
of SExtractor (Bertin & 



Arnouts 1996) as the representative galaxy magnitude 
Galactic foreground ext inction was co r rected using 

The 



(1998). 



the reddening maps from Schlegel et al. 
E(B — V) value applied to each cluster is 0.039, 0.042, 
0.028, and 0.016 for A119, A2670, A3330, and A389, re- 
spectively. Finally, fc-correction was performed by means 
of simple two-population (old plus young) modeling on 
the UV-optical continua from photometry consisting of 
GALEX FUV, NUV and CTIO it', g', r' bands. In 
the modeling, the old population has a fixed age of 12 
Gyr and the young component varies in age (0.01 Gyr 
< ty OU ng < 10 G yr) and mass fra ction (10~ 4 < f yoU ng < 
1). Section 3 in Yi et al. (20111 details the procedure. 
The average of the /c-correction values the for the cluster 
red-sequence galaxies are presented for each galaxy clus- 
ter in Table [3] The fc-corrected magnitudes are used in 
all of the analyses in this paper. 

The distance moduli used for the absolute magnitudes 
calculations are 36.25, 37.46, 37.88, and 38.32 for A119, 
A2670, A3330, and A389, respectively. 

3.2. Hydra spectroscopic data 

Hydra spectroscopic data were reduced using the 
hydra package in IRAF. First, the 2Kx2K images (af- 
ter binning along the dispersion axis) were pre-processed 
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Fig. 2. — Radial velocity histogram of the galaxy clusters. The 
Gaussian fits are plotted over the histograms. The dotted lines 
indicate the velocities of v r , center i 3cv r , with which cluster mem- 
bers are selected. The arrow in the panel (c) indicates a background 
cluster of Abell 3330 which appears in Section 4.2. 

to flat-fielding as with the imaging data. After the pre- 
processes, spectral apertures were traced from the two- 
dimensional images and fiber-to-fiber throughput varia- 
tion was corrected with a combined sky flat taken in each 
observing run. Wavelength calibration and sky subtrac- 
tion were also carried out. The procedures after pre- 
processing were performed via the dohydra task. All 
of the spectra were taken in three separate exposures 
for median combine. To combine the spectra, we nor- 
malized each spectrum using the continuum task, after 
which three spectra of a galaxy were combined into one 
spectrum with the scombine task. 

Galaxy radial velocities were measured from the com- 
bined spectra using the fxcor task in IRAF, which 
adopts a cross-correlation method between a rest -frame 
galaxy spectra l template and a galaxy spectrum ( |Tonry| 
fe; Davis||1979[ ). The early-type galaxy template and the 
luminous red galaxy template from SDSS were utilized 
for our radial velocity measurement , as we primarily tar- 
geted red-sequence galaxies in the Hydra observations. 

4. GALAXY SAMPLE SELECTION 

4.1. Cluster memberships 

Cluster memberships are determined by spectroscopic 
rcdshifts of galaxies. We measured the radial velocity 
dispersion (cr Vr ) of the cluster galaxies with the Hydra 
spectra. Because A2670 and the northern half of A119 
are observed by SDSS, we included the SDSS spectra of 
galaxies which we had missed in the Hydra observations. 
Radial velocity distributions from the combined catalogs 
(Hydra observations + SDSS archive) are shown in Fig- 
ure [2] The bin sizes are determined as one hundredth of 
the central velocities of the galaxy clusters. The center 
and dispersion of the radial velocity distributions are de- 
termined using a Gaussian fitting IDL code, gaussf it. 
The fitting result is plotted over the histogram with a 
solid curve. The central velocities and velocity disper- 
sions of the cluster galaxies are presented in Table |4j 
Cluster memberships are assigned to the galaxies in- 
cluded within ± 3<7„ r from the cluster central velocity. 
In Figure [2] the dotted lines mark the velocity ranges 
of the cluster memberships. As a result, we identified 
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Fig. 3. — Optical color-magnitude relationships of the four tar- 
get clusters with the MOSAIC 2 data. The dark-gray open circles 
indicate the spectroscopic members while the bright-gray open cir- 
cles present all of the galaxies from the field of deep optical im- 
ages. Red-sequences are defined by iterative linear fitting with 
the spectroscopic member galaxies. The dotted lines in each panel 
indicate the ±3(7 boundaries from the best fit. We defined the red- 
sequences with galaxies within the boundaries. The luminosity cut 
{M r i = —20) for this study is indicated with a gray dashed line in 
each CMR. 

133, 153, 130 and 114 member galaxies for A119, A2670, 
A3330 and A389, respectively. Among them, 32 and 39 
galaxies are supplemented from the SDSS spectroscopic 
data for A119 and A2670. 

To estimate -R200 fr om the velocity disp ersion, we 
adopted equation (8) of Carlberg et al. (19971. With the 
assumption of that density profile is a singular isother- 
mal sphere in the target clusters, -R200 can be defined 
as, 

The a is the characteristic velocity dispersion of a galaxy 
cluster. In our calculation, we employed the a Vr from the 
Gaussian fit of each galaxy cluster for the a. The Hub- 
ble constant, H n = 70km s _1 Mpc -1 , is applied to the 
equation instead of the Hubble parameter, H(z), assum- 
ing no evolution of H in the redshift range, < z < 0.1. 
As a result, -R200 is estimated as 1.53, 2.16, 1.84, and 
2.07 Mpc for A119, A2670, A3330, and A389, respec- 
tively. The MOSAIC 2 CCD does not cover the entire 
region of the galaxy clusters. The general width of a 
side of a stacked deep image is 36', which corresponds to 
1.9, 3.1, 3.6, and 4.4 Mpc at distances of A119, A2670, 
A3330, and A389, respectively. If it is assumed that a 
cluster center is located in the center of a deep image, 
the images cover 62% (0.95/1.53), 72% (1.55/2.16), 98% 
(1.8/1.84), and 106% (2.2/2.07) of i? 2 oo for each galaxy 
cluster. 

With R200, we estimated the mass of the clusters as 
well. All of the target clusters were more massive than 
2 x 10 14 M o . The mass of A2670 was - 7.3 x 10 14 M Q , 
proving that it is one of the most massive galaxy clusters 
in the universe. 

4.2. Red-sequence galaxies 

Because our interest is focused on the formation of 
massive red galaxies in galaxy clusters, we defined a red- 
sequence in each galaxy cluster first. The optical color- 
magnitude relationships of the target clusters are pre- 
sented in Figure [3j In the fi gure, the dark-gray open cir- 
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TABLE 4 

Spectroscopic properties of the clusters 



ID 


Vr, center 




z 


R200 


M200 


# of spectroscopic members 51 




(km/s) 


(km/s) 




(Mpc) 


(xlO 14 M ) 




A119 


1.32X10 4 


618 


0.044 


1.53 


2.6 


133 (32) 


A2670 


2.27X10 4 


871 


0.076 


2.16 


7.3 


153 (39) 


A3330 


2.67X10 4 


746 


0.089 


1.84 


4.5 


130 


A389 


3.37X10 4 


837 


0.112 


2.07 


6.4 


114 



a The numbers in parenthesis indicate the number of member galaxies identified using SDSS 
data. 
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Fig. 4. — The top panel shows the spectroscopic complete- 
ness of our data for the red-sequence galaxies along the r' band 
absolute magnitudes. Although the spectroscopic coverage de- 
creases as the galaxies become fainter, still more than 73% of the 
—21 < M r i < —20 photometric red-sequence galaxies are tar- 
geted. The bottom panel shows the spectroscopic member frac- 
tions among the photometric red-sequences along magnitude. The 
horizontal dotted lines mark 73% and 50% in the top and bottom 
panels, respectively, as the lowest fraction in each plot. With this 
plot, we learn that red-sequences are considerably contaminated 
by background/foreground galaxies even in the cluster fields. 

cles indicate the spectroscopic members while the bright- 
gray open circles present all of the galaxies from the 
field of deep optical images. Though the red-sequences 
are distinctively revealed in the CMRs, we empirically 
defined red-sequences using the spectroscopically iden- 
tified cluster members. However, it was not as simple 
as visual recognition of the sequences, as a red-sequence 
has a slight bend near the bright end of the sequence. 
Thus, we carried out two-step linear fitting by means 
of a 2- a clipping method. Initially, a preliminary fit 
was derived using data within a rectangle box which in- 
cludes galaxies near the bright end of an apparent red- 
sequence. Once the preliminary fit appears to be close 
to the actual sequence, we calculated the standard devia- 



tion (<7g 



of the data to the fit. Linear fitting was then 



performed again using the new data set confined by the 



new ±2cr„ 



color boundaries with 



g ,_ r > coior uouiiuaries wnn r < r bnghtest 
This two-step fitting process allowed us to obtain reason- 
able fitting results corresponding to what we observed 
in the CMRs. The zero points and slopes for the red- 
sequences in our targets clusters are 1.107 ± 0.079 and 
—0.026 ± 0.004, respectively. The small scatters are con- 
sistent with previo us findings of universal red-sequences 
in galaxy clusters (|Bower, Lucey fc Ellis||1992[ |Stanford^| 
Eisenhardt & Dickinson 199^T Using the final fits, the 
standard deviations were calculated again with the data 



within the ±2ct 



color boundaries. The final val- 



ues were oy_ r 
A2670, A3330 



•=0.028, 0.039, 0.036, and 0.037 for A119, 
and A389. The red-sequences are defined 



with galaxies within ±3ay_ r ' using the final fits and 
standard deviations. The color boundaries are expressed 
in dotted lines in Figure [3j 

We attempted to perform a complete spectroscopic sur- 
vey of the red-sequence galaxies during the Hydra ob- 
servations. The top panel in Figure [4] shows the spec- 
troscopic completeness of our data for the red-sequence 
galaxies along the r' band absolute magnitudes. Al- 
though the spectroscopic coverage decreases as the mag- 
nitude becomes fainter, still more than 73% of the —21 < 
M r > < —20 photometric red-sequence galaxies are tar- 
geted. The bottom panel shows the spectroscopic mem- 
ber fraction in the red-sequences along the magnitudes. 
With the plot, we learn that red-sequences are consid- 
erably contaminated by background/foreground galax- 
ies even in the cluster fields. In the case of A2670, 
though it has very high spectroscopic completeness for 
the faintest magnitude bin, the spectroscopic member 
fraction is in good agreement with the results from other 
clusters. Abell 3330 shows the low member fraction for 
the brightest magnitude bin in the bottom panel of Fig- 
ure 4] It is due to the brightest cluster galaxy of a close 
background cluster in the field. The background cluster 
is indicated with an arrow in Figure [2] (c) . 

The absolute magnitude cut of M r < —20 is applied for 
the galaxy sample selection because beyond that, spec- 
troscopic survey completeness is low and visual inspec- 
tion is difficult for the galaxies at z ~ 0.1. Finally, our 
galaxy samples are defined by the spectroscopic mem- 
bers in the red-sequences; they are also brighter than 
M r = —20 in each galaxy cluster. We term the galaxy 
samples 'RS sp ' throughout this paper. 

5. VISUAL INSPECTION 
5.1. The scheme 

Our approach to find post-merger signatures involves 
a visual inspection of the galaxies with the MOSAIC 2 
r' band deep images. The faint surface brightness level 
of fi r ' <~*j 30 mag/arcsec 2 is acquired in the deep images 
for all target clusters. The general Petrosian diameter 
along the major axis of the M r < = —20 galaxy is ~ 14.4" 
(~ 60 pixels) in the r' band deep images, corresponding 
to ~ 27 kpc at z = 0.1. Generously assuming that a 
typical seeing is 1 arcsec, the spatial resolution is ~ 1.9 
kpc at this distance. This can be regarded as the worst 
condition for our visual inspection. 

We attempted to establish the classification crite- 
ria for the visual inspection in a manner faithful to 
how the galaxies appear in the images. Therefore, the 
galaxies were classified as "Elliptical/SO" (E), "Post- 
merger" (PM), "Interacting" (I), "Pair" (P), "accompa- 
nying Faint Companions" (FC), and "Spiral" (S). "El- 
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Fig. 5. — This figure shows galaxy examples along the classification scheme. The first two columns show post-merger (PM) galaxy 
examples and the examples of interacting galaxies (I), and pair (P), and faint companions (FC) are presented in the following columns. 

liptical/SO" denotes the normal elliptical and lenticu- 
lar galaxies without any disturbed features. Addition- 
ally, "Spiral" indicates normal spiral galaxies. "Post- 
merger" denotes galaxies showing disturbed features, 
e.g., asymmetric structures, faint features, discontinuous 
halo structure, rings and dust lanes. These are consid- 
ered as galaxy merger remnants. "Interacting" galaxies 
also exhibit disturbed features as well as PM galaxies 
but they are close companions and interact with each 
other, i.e., they are on-going merger candidates. "Pair" 
refers to galaxies with a close companion without any 
disturbed features. We do not know whether they are 
actual pairs or simply a line-of sight effect of physically 
separated galaxies. Though we classified these objects, 
we do not include them as interacting systems in this 
paper. The classification "accompanying Faint Com- 
panions" was proposed after the preliminary visual in- 
spection of the overall populations of the galaxies. In 
this process, we found many galaxies accompanying faint 
dwarfish objects, mostly showing optical colors similar 
to those of the main galaxies. The small companions are 
usually extended along the radial direction to the center 
of the main galaxy or along the major axis of the galaxy. 
We classified them separately, suspecting that they are 
on-going minor merger systems or merger remnants cap- 
turing tidal dwarf galaxies that formed during the merger 
event of the past. Due to the ambiguity of identifying the 
systems, however, we did not include these galaxies in the 
analyses in this paper, too. Examples of the 'PM', T, 
'P', 'FC galaxies are presented in Figure [5] It should 
be noted that only PM and I samples are used in our 
analysis. 



TABLE 5 

Fractions of the post-merger, interacting, pair, and 
faint companion galaxies among rs sd 



ID 


PM 


I 


P 


FC 


A119 
A2670 
A3330 
A389 


27% (13/48) 
28% (25/89) 
22% (13/58) 
22% (17/78) 


0% (0/48) 
8% (7/89) 
5% (3/58) 
4% (3/78) 


10% (5/48) 
3% (3/89) 
9% (5/58) 
8% (6/78) 


10% (5/48) 
6% (5/89) 
3% (2/58) 
5% (4/78) 


among the RS 3p galaxies within 0.5-f?200 


A119 
A2670 
A3330 
A389 


26% (8/31) 
29% (18/63) 

13% (4/30) 
23% (10/43) 


0% (0/31) 
10% (6/63) 
3% (1/30) 
5% (2/43) 


10% (3/31) 
5% (3/63) 
7%(2/30) 

12% (5/43) 


16% (5/31) 
8% (5/63) 
0% (0/30) 
7% (3/43) 
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5.2. The result 

Visual inspection was performed for the four Abell 
clusters following the scheme introduced in the previ- 
ous subsection. In Table [5j we present the fractions of 
the post-merger (PM), interacting (I), pair (P) and faint 
companion (FC) galaxies among the KS sp of each galaxy 
cluster. We found that the ~ 25 ± 3% of the RS sp galax- 



FlG. 6. — The fractions of post-merger galaxies are plotted along 
the absolute magnitudes. The numbers beside the points indicate 
the count of the PM galaxies used to derive the fractions. Al- 
though there is a large spread according to the small values for 
the bright galaxies, the overall tendency is consistent between all 
four target clusters. This figure shows that post-merger features 
are more common as red-sequence galaxies become brighter in a 
cluster environment. 
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TABLE 6 

A CATALOG OF POST-MERGER GALAXIES. 



ID 
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Fig. 7. — Post-merger galaxy samples from four Abell clusters. The first column of each galaxy is a pseudo-color image generated with 
u'g'r' deep images. The second column is a surface brightness map of the galaxy, revealing faint features which do not appear in the color 
image. 



10 



Sheen et al. 




Fig. 8. — The BCGs in our cluster samples are presented with 
r' band deep images. The contours are drawn for fi r i = 19, 20, 
21, 22 mag/ arcsec 2 . The inset image in right bottom corner of 
each galaxy shows the central region indicated by a rectangle box 
in each image. A119, A2670 and A3330 showed faint structures in 
their halos. Moreover, A3330 and A389 revealed multiple cores in 
their centers. This implies that even the BCGs continued violent 
galaxy mergers until the recent epoch. 

ies exhibit post-merger signatures in the galaxy clusters 
at z < 0.1. The fractions of the other classes were usu- 
ally less than 10% for each class. If we consider post- 
mergers and interacting galaxies together, the fractions 
of the disturbed galaxies will be 27%, 36%, 28%, and 
26% respectively in A119, A2670, A3330 and A389. On 
average, the fraction of galaxies related to recent galaxy 
mergers is ~ 30 ± 4% in the clusters at z < 0.1. 

The FOVs of the MOSAIC 2 deep images do not 
cover the entire area within R2 00 for A119 and A2670, 
as demonstrated in Section |4.1| Therefore, we derived 
the fractions again with the RS"7 P galaxies located within 
0.5i?2oo of each galaxy cluster. The results are presented 
in the bottom half of Table [5] Although we applied 
more consistent areas of clusters, the result did not differ 
much from the previous one. The fraction of post-merger 
galaxies within 0.5i?2O0 was ~ 24 ± 4% on average. The 
value slightly lower than that (<~ 25 ± 3%) in the pre- 
vious sample is most likely a result of A3330 showing a 
distinctively low post-merger fraction within the half of 

-^200- 

While most of the fractions in a certain class show sim- 
ilar values between the clusters, it is interesting that no 
interacting galaxies were found among the bright red- 
sequence galaxies of A119. This may have been affected 
by the relatively small sky coverage of A119 (~ O.6.R200) 
by the MOSAIC 2 deep images. Because galaxy in- 
teractions would be more plausible at cluster outskirts, 
it is possible that we missed some galaxies out of the 
FOV. Another possible reason is that A119 is the small- 
est galaxy cluster among our target clusters and therefore 
includes fewer bright galaxies. The derived size and mass 
of the clusters are presented in Table [3] If we assume 
that the fraction of interacting galaxies within 0.5i?2oo 
of A119 may be similar to the value for A3330 (~ 3%), 
which is the second smallest galaxy cluster among the 
clusters, we expect only < 1 galaxies within the area. 
Therefore, there is a very low chance to find the bright 
interacting galaxies in the given FOV of the A119 data. 



The fractions of post-merger (PM) galaxies are plotted 
along the absolute magnitudes in Figure [6j The numbers 
beside the data points indicate the counts of PM galax- 
ies used to derive the fractions. Although there is a large 
spread according to the small values for the bright galax- 
ies, the overall tendency is consistent between all four 
target clusters. This figure shows that post-merger fea- 
tures are more common in bright red-sequence galaxies 
in a cluster environment. 

Samples of the selected post-merger galaxies are pre- 
sented in Figure[7j In this figure, the first column of each 
galaxy image set is a pseudo-color image generated with 
u'g'r' deep images. The second column is a r' band sur- 
face brightness map of the galaxy revealing faint features 
which do not appear in the pseudo-color images due to 
the luminous galactic centers. 

In fact, the BCGs in all four target clusters revealed 
extended faint features or multiple nuclei in our deep 
images. We present four BCGs separately in Figure [8] 
with their r' band deep image and the surface bright- 
ness contours. The inset images show the center of the 
BCGs indicated by the rectangle boxes. A119, A2670 
and A3330 BCGs showed faint structures in their halos. 
In addition to that, A3330 and A389 BCGs revealed mul- 
tiple cores in their centers. This implies that even the 
BCGs continued violent galaxy mergers until the recent 
epoch. 

We provide a catalog of the post-merger galaxies in 
Table [6] with their observational properties. 

6. MORPHOLOGICAL EXAMINATION 

6.1. Bulge-to-total ratios and residual images 

Galaxy bulge-to-total ratios (B/T) were measured for 
the spectroscopic members of the galaxy clusters. Ratios 
are derived from the radial surface brightness profiles of 
the galaxies, as measured by the ellipse task in IRAF. 
Because the cores of bright galaxies are usually saturated 
in the deep images, we merged the two profiles from the 
shallow image and the deep image for the central part 
and the outer part of a galaxy, respectively. The con- 
junction is made at the radius where both profiles show 
the same surface brightness. 

We found the best fits of the galaxy profiles to the 
composite galaxy model of the de Vaucouleurs profile 
(Sersic index, n = 4) and the exponential profile. The 
galaxy surface brightness model can be expressed as, 

= I e exp{c x ((R/R e ) 1/4 - 1)) + I d exp(-R/R d ). 

(5) 

Here, R e is the effective radius of the de Vaucouleurs pro- 
file and I e is the surface brightness at R e . The constant, 
c, is set as a free parameter to be determined. R d is the 
disk scale length and the I d is the surface brightness at 
R d . The fitting is pe rformed using the IDL routines in 
the MPFIT package ( |Markwardt]|2008[ ) . Usually such a 
fitting procedure is very sensitive to the initial values of 
the fitting parameters. To increase the success rate of 
the fitting, galaxy profiles are fitted with the de Vau- 
couleurs profile first. The preliminarily derived values 
of I e , R e and c are then used as an initial guess for the 
more complex main fitting with the equation ^ and all 
the parameters are determined again as well as I d and 
R d . The galaxy B/T is calculated using the following 
formula: 
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Fig. 9. — Galaxy radial profiles and the fitting results are pre- 
sented. The open rectangles indicate that the ellipsoid is fitted 
successfully for the annulus by the ellipse task. The best fit of 
the surface brightness profile is depicted using solid lines over the 
data points. The de Vaucouleurs profile and the exponential pro- 
file comprising the best fit are also presented, in dashed lines and 
dashed-dotted lines, respectively. This figure shows that the galaxy 
profiles are robustly measured to /i r r ~ 30 mag/arcsec. 



B/T = 



0-28R 2 d I d 



(6) 



We present examples of the galaxy profile fitting in 
Figure [9| The best fit of the surface brightness profile 
is drawn over the data points using the parameters de- 
rived from the fitting. The de Vaucouleurs profile and 
the exponential profile comprising the best fit are also 
presented. The figure shows that the galaxy profiles are 
robustly me asured to fj, r / ~ 30 mag/arcsec. 

Figure 10 shows the result of B/T calculations for the 
cluster members. In the upper panel of each graph set, 
galaxies are expressed in different colors along their B /T 
values in the optical CMR. The high values of B/T of 
bright red-sequence galaxies and the low values of blue 
galaxies are consistent with the general expectation. In 
the bottom panels, average B/T values for RS sp are plot- 
ted along the absolute magnitude. As the upper panels 
have already suggested, it was found that the average 
B/T becomes smaller as the galaxies become fainter. 

We also computed the averages for the post-merger 
galaxies separately and presented them in red filled cir- 
cles with the standard deviations. The results indicate 
that the average B /T values for post-merger galaxies are 
slightly higher than the average of the rest of the KS sp 
galaxies in the magnitude range. 

We classified galaxies with B/T > .4 as ellipticals and 
SOs following the criteria applied in Somerville fc Pri- 




"V* H&frvi 

■ *J 





m" 



i-o Fa 
0.8 r" 

0.6 L - 

0.4 : 

0.2 : 

0.0 : 



0.8 r 



2 - 

r 
-24 



Fig. 10. — The derived B/T values are presented for the spectro- 
scopic member galaxies brighter than M T i = — 20 in each galaxy 
cluster, upper: The B/T are expressed with color-coded filled cir- 
cles in the optical CMR of each galaxy cluster. The dashed lines 
indicate color ranges of the red-sequence galaxies. The graphs con- 
firm that bright red-sequence galaxies are mostly bulge-dominated 
systems with large B/T values while the blue galaxies show small 
B/T values uniformly, bottom: The averages B/T of post- merger 
galaxies among RS sp galaxies are plotted along the magnitude with 
the red-filled circles. The error bars show standard deviations of 
the mean values. Blue open circles are the average B/T of RS S p 
galaxies except the PM galaxies, presented for comparison. The 
number of galaxies counted for each average value is indicated be- 
side the data point. We demarcate B/T = 0.4 in the plots with 
the dotted lines, with which E/S0 galaxies are distinguished from 
spirals in this work. The graphs indicate that most of PM galaxies 
are bulge-dominated E/S0 galaxies. 
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Fig. 11. — Pseudo-color images, two-dimensional galaxy models and model-subtracted residual images of galaxies are presented from left 
to right in the same scale for each galaxy. The images of an early-type galaxy are presented in the top panels for comparison. The residual 
in the middle panel shows an example of discontinuous (in other words, structured) halo of an early-type galaxy while the bottom panel 
exhibits an obvious merger remnant with a long tidal tail and the spiral structures. 

galaxies demonstrate the robustness of our B/T calcula- 
tion for early-type galaxies. Pseudo-color images, two- 
dimensional galaxy models and model-subtracted resid- 
ual images are presented from left to right. An image of 
an early-type galaxy is presented in the top panels for 
comparison. Middle and bottom panels show examples 
of post- merger galaxies showing asymmetric, disturbed 
structures in their residual images. 



mack ( |1999 ). The results show that a large fraction of the 
featured galaxies are bulge-dominated systems as 77% 
(10/13), 72% (18/25), 77% (10/13) and 59% (10/17) for 
A119, A2670, A3330, and A389, respectively. This im- 
plies that possible contamination from disturbed spirals 
is not significant in our visual inspection. As a result, 
the fractions of post-merger galaxies among the bulge- 
dominated systems (B/T > 0.4) are about 40% (10/25), 
34% (18/53), 36% (10/28) and 48% (10/21) for A119, 
A2670, A3330, and A389, respectively. On average, it is 
38 ± 5% in our cluster sample. If we apply a more strict 
B/T criterion for bulge-dominated systems as B/T > 0.6, 
the average of post-merger galaxy fractions increases to 
46 ± 8%. 

Two-dimensional galaxy surface brightness models are 
constructed with the result of ellpse fitting and are sub- 
tracted from the deep images in order to check the fea- 
tures of the post-merger galaxies. The residual images of 
three KS S V ga laxies with similar B/T values are shown 



6.2. Asymmetry 

Using the r' band deep images, we measured the asym- 
metry index, A for the cluster members. The measure- 
ments were made within the area between 0.5Rp e t ro and 
2Rp etro of each galaxy avoiding galaxy centers usually 
saturated. 

The A can be expressed as, 



in Figure 11 At the same time, the images of three 



A = 



Eij\ni,3)-h8o(i,j)\ 



E,,-|i(*\i)l 



(7) 
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Fig. 12. — Asymmetry. A color composite image (first), surface brightness map (second), r band deep image (third), the same deep 
image rotated 180 degrees (fourth), and a normalized residual image (fifth). The green ellipses in the last column indicate 0.5i?p etro and 
2Rp etro of each galaxy in which A was measured. 

Prior to the calculations, all nearby objects around the 
galaxies were masked from the images as shown in the 
images with white masks. 

The average A of normal elliptical and spiral galaxies 
and the post-merger galaxies are plotted along their B/T 
values in Figure [13] While the figure confirms that late- 
type galaxies show large asymmetry values compared 
to early-type galaxies, we found that the asymmetry in 
the post-mergers arc slightly larger than that of normal 
early-type galaxies on average. Although the result cor- 
responds to simple expectation, however, the averages 
of A contain a large spread of the values as expressed 
with error bars in the figure. We suppose that the post- 
merger features in such a faint surface brightness level of 
our r' band data were not easily recognizable with the 
asymmetry index, A. 
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Fig. 13. — Asymmetry distribution of red-sequence member 
galaxies along their B/T. Compared to normal elliptical and spiral 
galaxies, PM galaxies showed a slightly larger A on average. 

where / is an original image and the Jiso is the image 
rotated by 180°. We demonstrate the calculation pro- 
cess with some examples of RS sp galaxies showing similar 
B/T values in Figure 12 The first column shows pseudo- 
color images and the second one displays r' band surface 
brightness map of the galaxies. The third and fourth im- 
ages show the original deep image and the same image 
rotated 180 degrees. The rotated images are subtracted 
from the original images and the normalize by the orig- 
inal images. The last column presents the normalized 
residual images with green ellipses indicating 0.5Rp e t ro 
and 2Rp etro of each galaxy in which the A was measured. 



7. DISCUSSION 

We investigated the post-merger signatures of red- 
sequence galaxies in rich Abell clusters at z < 0.1. Deep 
images in u', g' , r' and medium-resolution galaxy spec- 
tra were taken for A119, A2670, A3330 and A389 with 
a MOSAIC 2 CCD and a Hydra MOS mounted on the 
Blanco 4-m telescope at CTIO. Post-merger signatures 
were identified by visual inspection of their disturbed 
features, e.g., asymmetric structures, faint features, dis- 
continuous halo structures, rings and dust lanes. Most 
(~ 71%) of the featured galaxies were found to be bulge- 
dominated systems (E/S0 galaxies) from the radial sur- 
face brightness profile fitting. On average, the asymme- 
try of the post-merger galaxies is always slightly larger 
than the mean of normal galaxies within the same bulge- 
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TABLE 7 

Comparisons between cluster and field 





Class 


Cluster 


Field a 




PM 


25 ± 3% 


35% 


Red b 


I 


5 ± 1% 


18% 




Total 


30 ± 4% 


53% 




PM 


38 ± 5% 


49% 


Bulge-dominated 


I 


4 ± 1% 


21% 




Total 


42 ± 6% 


70% 



a The fractio ns for the field envi ronment were 
adopted from Ivan Dokkum|p005t . 
b Fractions for the cluster are derived with the 
RS S p galaxies in this paper, while the field red 
galaxies are denoted with B — R colors. 
c Fractions with only bulge-dominated galaxies 
among the red galaxies. For the cluster, galax- 
ies with B/T > 0.4 are included while visually 
classified F/SO galaxies are considered for the 
field. 

to-total ratio range. 

We noted that cluster red-sequence galaxies went 
through galaxy mergers until the recent epoch. In this 
work, — 25% of the bright (M r < -20) cluster red- 
sequence (RS sp ) galaxies exhibited post-merger signa- 
tures in four target clusters consistently. The fraction in- 
creases to ~ 38% with the bulge-dominated RS sp galaxies 
(B/T > 0.4). In addition, brighter galaxy samples show 
higher fractions of post-merger signatures in the clusters. 
In Table 7j we compared the average of fractions from 
bell clusters with the result of field environ- 
The clusters and field 



the four A 

ment from van Dokkum ( 2005 1 
environments were compared using two subsets of galaxy 
samples, red galaxies and bulge-dominated red galax- 
ies. Our RS sp galaxies were compared with the l umi- 
nous red galax ies selected using B — R colors in |van| 
Dokkum (2005), and the bulge-dominated galaxies (B/T 
> 0.4) were compared with the visually classified field 
E/S0 galaxies in the same paper. We presented the 
fractions of post-merger galaxies and interacting galax- 
ies separately as well as the combined fractions which 
indicate the fraction of disturbed galaxies related to past 
or ongoing galaxy mergers. The van Dokkum sample 
does come with spectroscopic rcdshift; spectroscopic red- 
shifts were available only on 5 galaxies out of 126 and so 
no k-correction was applied. Based on a small subsam- 
ple he argues that his galaxies are roughly at redshift 
of 0.1, which is comparable to ours. The colors of his 
bulge-dominated galaxies are also comparable to ours af- 
ter proper color conversions using Lupton (2005). 

The table demonstrates that the overall tractions in 
clusters are lower than those in the field. We find it dif- 
ficult to understand the "interacting" galaxies in these 
clusters. They appear to be beginning their interac- 
tion/merger, and so we classified them as "interacting" 
or pre-merger. Considering that interacting systems with 
small companions are likely missed in our search, it is 
likely that our fraction for interacting galaxies is a hard 
lower limit. For example, some of our "faint companion 
galaxies" (which were excluded in our analysis) could be 
"interacting". Direct satellite-satellite mergers are sup- 
posed to be extremely rare in large halo environments, 
and so it is difficult to understand them. 

The fractions of post-merger galaxies are not very dif- 
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Fig. 14. — Fractions of PM galaxies along the projected distances 
from the BCGs. The size of distance bin is 0.2i?2oo m each galaxy 
cluster. Contrary to the general assumption, which would predict 
more frequent mergers at the outskirts of a cluster, the fraction 
of PM galaxies does not change much along the clustocentric dis- 
tances. 

ferent between the cluster and the field. This is not 
compatible with a simple theoretical expectation based 
on merger timescale as a function of relative speeds of 
galaxies. Therefore, we suggest that many massive red- 
sequence galaxies formed in a less dense environment 
through galaxy mergers and entered into the cluster 
through dark matter halo mergers. 

An important caveat exists on this comparison between 
van Dokkum's work and ours. Merger timescales and 
thus probability heavily depends on mass ratio of collid- 
ing systems, but it is very difficult to extract the mass 
ratio information from observed images. We thought it 
would be the best attempt to use similar instruments for 
similar depth (exposure) on galaxies at similar distances. 
But even in this case, we are not free from a possible en- 
vironmental (cluster vs field) dependence on mass ratios 
of galaxy mergers. 

The radial distribution of the post-merger galaxies in 
the clusters provides another hint supporting this sce- 
nario. The fractions of post-merger galaxies are derived 
along the projected distance from the BCG in the scale 
of i?200 of each galaxy cluster, as shown in Figure [l4j 
We found that the fractions do not change much along 
the clustocentric distances. In situ mergers can scarcely 
explain this result because it is difficult to expect fre- 
quent mergers in the central region of clusters compared 
to the outskirts. This implies again that a large fraction 
of massive red-sequence galaxies with post-merger signa- 
tures likely merged in a less dense environment and then 
fell into the central region of the cluster. 

Is it possible for a galaxy to maintain the post-merger 
features produced in the outskirts until the galaxy arrives 
at the cluster center? It has been suggested for a long 
time that tidally ejected mat erials during galaxy mergers 



may return over many Gyr (Hernq uist fc Spergel 1992 
Hibbard fc Mihos|1995[ ). |Lotz et al.| ( |2008p also proposec 
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Fig. 15. — Merging timescale of galaxies in a cluster of 5xlO 14 M 
of our SAM along the mass ratio. If we assume that a central galaxy 
is the BCG, bright satellite galaxies (M sa t > 0.25M cen ) merge into 
the BCG within ~ 4 Gyr at z < 0.5. 

that the features of interacting galaxies can be observable 
for up to ~ 2 Gyr by applying stellar population models 
to the simulated equal-mass gas-rich disc mergers. We 
note that the duration of merger feature is highly sen- 
sitive to galaxy type, merger geometry mass ratio, and 
last but not least the observing surface brightness limitQ 
The dynamical friction timescale of a satellite galaxy 
to the central cluster halo would be another parameter 
to compare with the observable timescale of post-merger 
features. In this work the dynamical friction timescale 
is defined as the timescale between the epochs of when 
the halo of satellite galaxy is started to merge into the 
cluster halo at the cluster virial radius and when the 
satellite galaxy have finally merged into the BCG in the 
model. If a dynamical friction timescale for a galaxy 
is much longer than the timescale of the post-merger 
features, our results should be interpreted in a differ- 
ent light. By taking advantage of a semi-analytic ap- 
proach, we obtained the relationship between the merger 
mass ratio and the merger timescale of galaxies to the 
BCG. Our semi-analytic model (SAM) is based on N- 
body backbone dark matter merger trees and physical 
ingredients that govern the baryonic properties of galax- 
ies. The merger times cales in the model a re calculated 
using the formulae of Jiang et al. (2008), which pro- 
vides a modified form of Chandrasekhar's formula. Given 
that Chandrasekha r's formula generall y underestimates 



merger timescales, Jiang et al 
into merger timescales derived 



([2008]) fits the formula 
rom hydro dyna mic sim- 



ulations (see also Kimm, Yi, & Khochfar 2011 ). 

We selected a cluster in the model which has a clus- 
ter mass, 5xl0 14 M Q , comparable to our cluster sam- 
ples. Figure [15] shows the merging timescales of all satel- 
lite galaxy mergers along the mass ratio between the 
BCG and the satellite, consequently forming a cluster 
of 5xlO 14 M0 at present. The redshift of each merger 
was expressed using different symbols, showing that the 
merger timescale was shorter at a higher redshift for the 
same mass ratio in the model. This arises because the 
size of the cluster was smaller at an early epoch. We can 



1 The suggestion of 2 Gyr visibility of merger features by |Lotz et] 
|aL] l |2008| was based on the surface brightness, fi = 25 assumption. 
Trie visibility time is likely longer in our case because our imaging 
was deeper. 



assume that the oldest observable post-merger features 
in our target clusters are produced at z ~ 0.5, where the 
difference in the lookback time from z ~ 0.1 is about 
4 Gyr. The figure shows that the dynamical friction 
timescale of a satellite galaxy as massive as M sat /M cen > 
0.25 will be < 4 Gyr below the redshift, z = 0.5. If we 
assume that the magnitude difference between the BCG 
and the second BCG is 1 magnitude, then the mass ratio 
would be ~ 0.4 (M sat :M ce „ = 1:2.5). The figure indi- 
cates that then the second BCG will merge into the BCG 
within 3 Gyr at z < 0.5. Therefore, it appears that mas- 
sive post-merger galaxies can fall into the cluster center 
and maintain their disturbed features, even in low red- 
shift (z < 0.5). The consistent fractions of the featured 
galaxies along the clustocentric distance most likely arise 
from on-going halo mergers which started during various 
epochs. 

The average count of post-merger galaxies in our clus- 
ter samples was 17 ± 6 galaxies (13, 25, 13 and 17 galax- 
ies for A119, A2670, A3330 and A389, respectively). We 
also counted the number of halo mergers related to the 
current cluster halo in the SAM, to check whether a clus- 
ter can experience sufficient halo mergers to absorb mas- 
sive galaxies since z = 0.5. In the model of a 5xl0 14 M Q 
cluster, we found 21 satellite halos merging into the clus- 
ter at z < 0.5, with galaxies brighter than M r > = —20. 
It theoretically supports that a massive cluster may have 
enough halo mergers to take recently merged galaxies 
in at recent epoch. In this preliminary analysis, we as- 
sumed that merger features from the previous halo en- 
vironment last roughly the same time as the subhalo's 
dynamical friction timescale. This is certainly an over- 
simplification. A more realistic calculation would require 
accurate merging halo mass ratios, merging galaxy mass 
ratios, and merger feature timescales throughout cluster 
evolution history. 

In conclusion, we found that ~ 42% of massive, bulge- 
dominated red-sequence galaxies in galaxy clusters have 
continued their mass assembly through galaxy mergers 
until the recent epoch, as ~ 70% of the field galaxies have 
done. Although the fractions of post-merger galaxies in 
clusters are lower than that in the field, it is still too high 
compared to the expectation from apparent fractions of 
galaxy interaction in galaxy clusters. Therefore, we pro- 
pose that most of those post-merger galaxies were as- 
sembled in a low-density region and fell into the current 
cluster via halo mergers. We have supported this sce- 
nario with theoretical predictions using a semi-analytical 
model. 

We do not speculate the progenitors of the post-merger 
galaxies in this paper, regardless of whether they are 
remnants of dry mergers or wet mergers. This will be 
investigated in depth in an upcoming paper which will 
discuss on the ultraviolet properties of the galaxies most 
likely related to the residual star formation induced by 
galaxy mergers. 
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